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Theory of a weakly non-ideal Bose gas in the canonical ensemble is developed without assumption 
q ' of the C-number representation of the creation and annihilation operators with zero momentum. 

' It is shown that the pole of the "density-density" Green function exactly coincides with the 

, Bogolybov's phonon-roton spectrum of excitations. At the same time in the one-particle excitation 

spectrum a gap exists. The value of this gap is connected with the density of particles in the 
Z\ ' "condensate". 

Q 

oo 



in 



X 



PACS number(s): 05.30.Jp, 03.75.Kk, 03.75.Nt, 05.70.Fh 



I. INTRODUCTION 

1. Starting with the Bogolyubov's papers [1,2], the microscopic theory of the degenerate Bose gas has been based 
on the special assumption that the creation and annihilation do operators of particles with zero momentum can 



be replaced by a C-number 

4 = a = (< N >) 1 / 2 , (1) 



where < N > is the average quantity of particles in the state with p — ("condensate"). 

This assumption leads to necessity of introducing in the theory the anomalous averages ("quasi-averages"), which 
^ | is unnatural for the homogeneous and isotropic system under consideration. Though the Bogolyubov theory (see, 
*73 ■ e.g., [3,4]) gives rise to many important and widely recognized results, such as the expression for the spectrum of 
excitations, explanation of the experimental data and agreement with the Landau superfluidity condition, there are 
serious doubts on the validity of the relation as well as on the agreement of the Hamiltonian corresponding to the 
assumption |T]) with the original Hamiltonian of the system under consideration [5,6]. 

In this work we suggest a selfconsistent noncontradictory description of the degenerate Bose gas in which we do not 
use the assumption (!]). In this way we show, in particular, that the gap in the spectrum of the one-particle excitations 
exists. The statement that the spectrum of one-particle excitations has a gap in parallel with the usual phonon-roton 
branch of the collective excitations in degenerate Bose gas has been formulated in [7] . Later the existence of a gap has 
been suggested in [8]. Spectra, thermodynamics and dynamical structure factor for weakly non-ideal Bose gas below 
ly-^ . the condensation temperature T < Tq have been considered in [9-11] in terms of the dielectric formalism, similar to one 
in plasma-like systems. The dielectric formalism on the basis of the Bogolyubov assumption |T]) for the operators ajj" 
and do was developed in [12]. Recently the gap-existence problem for a weakly non-ideal Bose gas has been considered 
by different methods in papers [13-16]. In fact, the paper [14] confirmed (independently) and developed not only the 
| statement of Ref. [7] about existence of a one-particle gap in the weakly non-ideal Bose gas below the condensation 
temperature, but also the results of the papers [9-11] on a possibility of considering the superfluid system without 
symmetry breaking. Of course, the existence of a gap disagrees with a standard opinion that the gap is missing [3]. 
But we do believe that the conclusion about a missing gap is related to the Bogolubov simplifying assumption (jTJ, 
$— i ' and is not obliged to persist in a more rigorous theory. 

It is necessary to stress that in the early papers both Landau and Bogolyubov admitted that the one-particle 
excitation spectrum can have a gap, but later they dropped this idea, because of contradiction with observations of 
the phonon-roton branch of spectrum at small q in the neutron scattering experiments. The approach which includes 
the both spectra, with and without a gap, has not been found and coexistence of these two branches of excitations 
was not suggested. 

We suggest that in the limit of strong degeneration T — > 0, where T is the temperature of the system, all the 
particles tend to occupy the zero-momentum state, which means that 

lim < N a >=< N >, (2) 



where < N > is the average number of particles in the system. This statement is confirmed in the paper by the 
self-consistent consideration. Below we use the canonical ensemble where N is a given C-number. 
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II. STATISTICAL SUM AND AVERAGES IN CANONICAL ENSEMBLE 

Let us consider the Hamiltonian of a non-ideal Bose gas of particles with zero spin and a mass to in a volume V 

H = £ P a t a P + 7^7 U (9) a p 1 -q/2 a p 2 +q/2 a P2-q/2 a Pl+q/2 5 (3) 

v ~" q,pi,P2 

where a+ and a p are the creation and annihilation operators of particles with the momentum hp, 

[a P2 ,ap l ] = a P2 a p " l — ap" l a P2 = <5p!,p 2 , (4) 

e p = h 2 p 2 /2m is the energetic spectrum of a free particle and u(q) is the Fourier-component of the inter-particle 
interaction potential. 

It is convenient to extract the particular term Uq with q = from the sum over q in the Hamiltonian ((3j). This 
term can be written as 

^ 0=U ( 0) 2K ' N = H a P a t» «(<>) = «(« = 0) = u(g 0), (5) 

p 

where TV is the operator of total number of particles, and N is a C-number N =< TV >. Here and below the brackets 
< ... > mean the canonical-ensemble averaging. Since Uq is also the C-number one can write the statistical sum of 
the system under consideration in the form 

Z = Sp exp(-H/T) = cxp (-11(0) ^^" - ) Sp exp(-H /T), (6) 
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where 



H = H - U = £ p «p a p + ^? zZ u (9) a p 1 - q /2 a p 2 + q /2 a P2-q/2 a Pi+q/2, (7) 
p q#o,pi,p 2 

Therefore, to provide convergence of the statistical sum J7]) in the thermodynamic limit V — > 00 N — > oo, n = N/V = 
const., the known condition u(0) > has to be fulfilled. As Uq is a C-number, it does not affect any averaging at all, 
and for calculation of average values the Hamiltonian H is equivalent to H . For an arbitrary operator < A > we get 

< A >= Z- 1 5p{exp(-i?/T)vl} = Z^SpieM-Ho/^A} =< A > ; Z Q = Sp cxp(-H /T) (8) 

Let us now consider a more complicated situation, connected with calculation of a time dependent correlation 
functions fit) of the type 

f(t) =< [A(t),B(0)} >, A{t) = exp(iHt/h)Aexp(-iHt/K). (9) 
The time dependence of the operators a+(£) and a p (t) can be represented in the form 

a+(t) = exp{iH t/h)a+ exp(-iH t/H) exp(iNu(0)t/Vh). (10) 

a p (t) = cxp(-iNu(0)t/Vh)exp(iH o t/h)a p exp(-iH o t/h). (11) 

Therefore, if in each of the operators A and B the amounts of creation and annihilation operators coincide (which is 
typical for the operators of the physical variables) the time-dependent correlation function can be written as 

f(t) =< [A{t),B(0)] > , A(t) = exp(iH t/h)Aexp(-iH t/h). (12) 

On this basis below, in the framework of the canonical ensemble, we consider the average values with the Hamiltonian 
Hq ([7]). The free energy of the initial system with the Hamiltonian H (|3j), according to 10- (7]) reads 

F = -TlnZ = U + F , F = -TlnZo- (13) 
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III. EQUATIONS FOR "DENSITY-DENSITY" GREEN FUNCTION 

Experimentally the spectrum of collective excitations is found usually from data on the well observable maxima in 
the dynamical structure factor 5(q, u) for 0, 

1 f°° 

S(q,u) = — exp(wt) < p q (t)p_ q (0) > dt, (14) 

* J — OO 

PqW = a p-q/2W a p+q/2<X)> ( 15 ) 

P 

where p q (<) is the Fourier-component of the operator of particle density in the Heisenberg representation. The 
dynamical structure factor S(q,u>) (fl4"l) is directly connected [17] with the retarded density-density Green function 
X H (q, z) which is analytical in the upper semi-plane of the complex variable z (Imz > 0), 

S(<1^) = —, f h /Tn WW + »0), (16) 

1 — cxp(— nu/1 ) 

i f°° 1 
X fl (q, z) = -— dt exp(izt) < [p q (t)p_ q (0)] > = — << p q | p-q »„, (17) 



The definitions (fl6|) . fT7j) have to be taken in the thermodynamic limit, where 

lim < iV > =< N > = N. (18) 

According to Eq. (fT4|) the function x(?i z ) can be represented in the form 

x(q, z) = — Y f (p, q> z ), f (p, q> z ) = « a p- q/2 a P +q/2 1 p-q))^ ( 19 ) 
p 

The equation of motion for the function F(p, q, z) with the Hamiltonian Ho, determined by (Eq. Q), can be written 
in the form 

(kz + £p_ q /2 - ^p+q/2) F (P> <l> Z ) = /p-q/2 ~ /p+q/2 - 

y H U ( k )J2(( ( a p+k-q/2 a p 1 -k/2 a Pi+k/2 a P+q/2 - a p- q/2 a p 1 -k/2 a Pi+k/2a P -k+q/2) I P-q))z (20) 
k pi 

Here / p is the one-particle distribution on the momenta hp function 

/ p =< a+ap > (21) 

For a temperature T < Tq, where To is the temperature of condensation the one-particle distribution function / p can 
be represented as [18] 

f P =< N a > S pfi + /J(l - S Pi0 ), (22) 

where -/V = a^ao is the operator of the quantity of particles with the momentum equal zero ("condensate") , 
fp =< a+a p >o is the one-particle distribution function with non-zero momenta (the "overcondensate" states). 
Therefore, 



™ = «o + ^£/p T = ^+ / t^/J, (23) 



where no =< No > /V is the average density of particles in the condensate. From ([20)1 and (|22|) we find that 
the function F(p, q, z) has singularities at p = ±q/2. Therefore, the density-density function x(g, z) (fT9|) can be 
represented in the form 

X(q,z) = ~F(q/2,< lt z) + ^ F(-q/2, q, z) +1 ]T F T (p,q,z) (24) 
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The index T means, that the respective function describes the "overcondensate" particles. The singularities 
(conditioned by the condensate) in this function are absent. Then in the last term in Eq. lf24|) we can change summation 
by integration over momenta. The functions F(±q/2,q, z), extracted above, satisfy, according to Eqs. ((20|) . (|22|) . to 
the exact equations of motion 

(hz ~ e q ) F(q/2, q, z) = [(N Q ) ~fq]-y^2 u i k ) («k «p 1 -k/2 a Pi+k/2a q - 

fc#0 Pi 

a (t a p 1 -k/2 a Pi+k/2 a q-k) I P-q))z (25) 

(hz + e q ) F(-q/2,q,z) = -[(N ) - /J] - — ^ u(k) (a^a+^^a^+^ao - 

fc#0 pi 

a -q a p 1 -k/2 a Pl+k/2«-k) | P-q))z- (26) 

Let us consider further the case of strongly degenerate gas, where T — > 0. To find for this case the main terms, following 
to the Bogolyubov's procedure, let us extract in the right sides of Eqs. (|2o]) . f[26|) the terms, which are determined by 
the maximum quantity of the operators and ao. In the limit of strong degeneration, taking into account ( fl8|) . we 
can omit the other terms in the course of calculation of F. Then, since 

/q = /-q (q^ 0) Eqs. {M), ([Ml take the form 
(hz T e q ) F ( -°\±q/2, q, z) = ±[{N ) - / q T ] ± yu(q)(( (a+a+a^ao + a+al q a ao) I P- q ))z (27) 

Eqs. lf27f are exact at low temperature, since they determine the functions F^ ' (±q/2, q, z). For calculation the Green 
functions in the right side of Eqs. f27j) it is necessary to make some approximations. 



IV. DETERMINATION OF THE COLLECTIVE EXCITATIONS 

It should be emphasized that the Bogolyubov's approach, in which the operators aj and ao are considered as 
C-numbers leads to violation of the exact relations (|25|) , f26|) . In this approach {TJ , the main term F^ ) of the function 
F has the form [19] 

f (0) (±q/2, q, z) = (N Q ) (( a ±q | a+ q )) z (28) 

Respectively, instead the exact equations for two-particle Green functions ([25)) , l(26|) we obtain the equation of motion 
for one-particle Green functions form (|28"ll . As is shown below, the equations of motion for the two-particle Green 
functions without the approximation on the C-number representation of the operators n ao are essentially different 
from the equations for the one-particle distribution functions. Therefore, for calculation of the Green functions in the 
right part of l(27j) . the assumption {J} cannot be used. At the same time the idea on C-number approximation for some 
operators is, itself, very attractive. Below we use this idea in the variant alternative to the Bogolyubov assumption. 
According to fl8|) it is natural to accept that for calculation of the Green functions in the limit of strong degeneration 
T — > the C-number approximation has to be applied not to the operators a q h ao, but to the operator of the number 
of particles in the "condensate" Ao 



N = (N ). (29) 



In this case from lj27f directly follows 



(hz - e q ) F<°> (q/2, q, z) = [(No) - ft] + (q/2, q, z) + F^ (-q/2, q, z)} (30) 

(hz + e q ) F<°) (-q/2, q, z) = - [(No) - ft] - ^^u(q){F^ (q/2, q, z) + F^ (-q/2, q, z)} (31) 
From Eqs. i[30|) . (|3ip one can find the solutions for the functions F(±q/2, q, z) 

m /9 _ [<N >-fT](hz + e q ) [<N Q >-fZ](hz-e q ) 

F(q/2 ' q ' nz) = (hzy - (hu( q )Y ; *W 2 .**) ( hz f (hu( q )Y (32) 
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fiu(q) = ^e 2 q + 2n u(q)e q (33) 

The relation ([33j) for the spectrum hw(q) corresponds exactly to the known Bogolyubov expression [1,2]. By 
substituting to (fl"9"[) , and taking into account that for the case of strong degeneration the contribution of the functions 
F T (p, q, z) is negligible, we obtain the expression for the main term X (Q> z ) of the "density-density" Green function 

As it is well known, the singularities of the function x R (q> z ) determine the spectrum of collective excitations in the 
system. Therefore, under the assumption about C-number behavior of the operator No we obtain the Bogolyubov's 
result for the spectrum of the collective excitations in the degenerate and weakly interacting Bose gas. However, the 
question on the term / T / < No > in the figure brackets of ([33[) still exists. The problem is in the behavior of the 
function for the ideal Bose gas [20] 

/^={ex P (^)-ip (35) 

In the limit of small wave vectors q the function f^ 1 converges at non-zero temperatures (further in the text 1/q 2 - 
divergence). Moreover 



lim lim C ^ lim lim ft a (36) 



The similar problem arises when (|22]) is used. 



V. ONE-PARTICLE EXCITATIONS, GAP AND THE SELF-CONSISTENT DISTRIBUTION 

FUNCTION 

To calculate the distribution function f q for the "overcondensate" particles in Bose gas let us consider the one- 
particle Green function g R (q, z) 

g R (q 1 z) = ((a cl \a+)) z , q/0 (37) 



This function is directly connected with the distribution function by the relation [21] 

—g<{q,u}l g < (q,w) = -2hlexp( — J - 1 j Im g R (q,LJ + .0) (38) 
The equation of motion for the Green function g R (q, z) for q ^ reads 

(hz-e q ) g R (q,z) - l + -^^w(fc)^((a+ +k a P a q+k | a+)) z . (39) 

fe^O p 

As for the "density-density" Green function, we consider the case of a strong degeneration and extract in the right- 
hand side of Eq. ([39]) the terms with maximum quantity of the operators a q and ao . Then from Eq. ([391 we find 

(hz-eq) g R (q,z) = l + —u(p){(a^a <l a \ a+)) z (l - 5 q , ), (40) 
Further, by using again the assumption about the C-number behavior of the operator No, from Eq. (|40[l we obtain 

^ = fe^V (41) 

The expression for the spectrum of the one-particle excitations E q is given by 

E q = e q + n u(q), (42) 



From Eqs. (071) , Q4QJ , gl]) for T < T we obtain 

t 1 

f " = eME g /T)-r (43) 
Therefore, the function /J is finite for q — > 0. Moreover, in the limit of strong degeneration T — > 

/ 9 T - (44) 

for arbitrary values of g, in contrast with the case of of the ideal Bose gas. Therefore, the representation lj22|) for the 
one-particle distribution function f p is valid, the initial suggestions ([2| . (fl8|) are satisfied and the used procedure of 
extraction of the main terms is self-consistent. 

As it follows from l!42|) in the spectrum of one-particle excitations the appeared gap equals to 

A = E q ^ = no«(0). (45) 

A value of the gap is determined by the density of particles in the "condensate". The existence of the gap permits 
to extend essentially the applicability of the results obtained for T — ► 0. It is obvious, that in many applications the 
condition T — > is equivalent to the condition T -C A. 

Taking into account |44|) we can rewrite Eq. l[34"ll for the function x (°j z ) m the form 

- VB^fas? (46) 

On the basis of 1|46|1 practically all known results for the thermodynamical functions of the degenerate Bose gas can 
be reproduced, as it was done in [9-11]. 

Therefore, in contrast with the approach, based on the C-number approximation for the operators Oq and ao, by 
use the C-number approximation for the operator N we find the spectra of the collective and one-particle excitations 
are different. Both spectra, as it's easy to see, satisfy the Landau condition for superfluidity. For the one-particle 
spectrum the Landau condition is satisfied for the transitions between the "condensate" and the "overcondensate" 
state. The Landau condition is, naturally, violated for transitions between the "overcondensate" states. 



VI. CONCLUSIONS 



Summarizing the performed consideration we can assert that by calculation of the Green functions for the highly 
degenerate Bose gas on the basis of the C-number approximation for the operator Nq : 

A) The problem of 1/g 2 divergence, which arises for the ideal Bose gas, can be solved; 

B) The system has two different branches of excitations - the one-particle and the collective ones, both satisfying 
the Landau condition of superfluidity; 

C) The one-particle spectrum of excitations contains the gap in the region of small wave vectors, connected with 
the existence of the "condensate" ; 

D) The spectrum of the collective excitations corresponds to the "phonon-roton" excitations observed in the 
experiments on the inelastic neutron scattering; 

E) The necessity of anomalous averages (quasi-averages) for the description of the degenerate Bose- gas is absent. 
Therefore, the inter-particle interaction in Bose-systems leads not only to the drastic difference (in comparison 

with the ideal Bose gas) in the structure of the collective excitations, which are described by "density-density" Green 
function, but also to the crucial change in the distribution function of the one-particle excitations and in the one- 
particle spectrum of excitations for the "overcodensate" particles. 

On the basis of the obtained results the special diagram technique can be developed, similar to [19], but with the 
use of the C-number approximation for the operator Nq. 

The principle difference between the results of this work and the results of the "traditional" C-number approximation 
for the operators and ao is the existence of the gap in the spectrum of the one-particle excitations. The above 
analysis shows that this gap cannot manifest itself in the experiments on nonelastic neutron scattering in superfluid 
Helium [22,23]. However, such possibility cannot be excluded [16] in the experiments on the Raman light scattering. 
Moreover, in [24], where such experiments are described, there is a direct indication of existence of the gap. 
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